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Abstract. The class of nonlinear evolution equations (NLEE) - gauge equivalent to the N-wa.ve equations 
related to the simple Lie algebra g are derived and analyzed. They are written in terms of S{x,t) £ g 
satisfying r = rankg nonlinear constraints. The corresponding Lax pairs and the time evolution of the 
scattering data are found. The Zakharov-Shabat dressing method is appropriately modified to construct 
their soliton solutions. 
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1 Introduction and Preliminaries 

It is well known that the Lax representation (|l]) 

[i(A),Af(A)] =0 (1) 

is invariant under the group of gauge transformations. 

One of the first nontrivial examples of gauge equivalent 
nonlinear evolution equations (NLEE) is provided by the 
nonlinear Schrodinger equation (NLSE) |P,|l],||]: 

iut +Uxx + '2\u'^\u{x,t) = 0, (2) 

and the Heisenberg ferromagnet equation (HFE): 



^(°)(:r,i)=5<"^-W"(^,t); 



(3) 



obviously (5*^°^)^ = 1. The equivalence between (^ and 
(0) is based on the fact that g''^^{x^t) is determined by 
u(x,t) through (see |^): 

^^!^l + ,(o)(^,t)5(o)(^,i) = 0, (4) 



g(")(x,i) = 



dx 

u 

-u* 



lim g'^^^x^t) = 1. (5) 



Both equations are infinite dimensional completely in- 
tegrable Hamiltonian systems. The phase space A^nlse 
is the linear space of all off-diagonal matrices g^°^ (cc, t) 
tending fast enough to zero for x — *■ ±c». A hierarchy of 
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pair-wise compatible symplectic structures on A^nlse is 
provided by the 2-forms: 



(fc) _ 

NLSE ^ 



dxtr (6q'^°^ AA^[(73,Sq'^°\x,t)]] . (6) 



The phase space A^hfe of the HFE is the manifold of all 
S'*^°-'(a;, t) determined by the second relation in (^). The 
family of compatible 2-forms is: 



"hfe ^ A 



dxtr 



Ai^-[5W,,55W(a;,i)]) . (7) 



By A and A we have denoted the recursion operator of the 
NLS type equations and its gauge equivalent [Q . The spec- 
tral theory of these two operators underlie all the funda- 
mental properties of these two classes of gauge equivalent 
NLEE, for details see |4|. Note that the gauge transfor- 
mation relates nontrivially the symplectic structures, i.e. 



{k+2) 



NLSE — '^HFE d'Sl- 

The NLSE is solvable by the inverse scattering method 
applied to the Zakharov-Shabat system. It can be gener- 
alized to any simple Lie algebra g of rank r > 1 by 

""To|,|Til: 



(8) 



L{\)i; = + q{x, t) - XJ ] tpix, t, A) = 0, 



where q{xA) and J are elements g. A natural choice of the 
gauge in (||) consists in choosing J to be a constant regular 
element of the Cartan subalgebra f) of g. Next by a gauge 
transformation L{\) — > g^^L{\)gQ where go{x,t) G <3 
commutes with J we can eliminate the diagonal elements 
of ^(a:, t) and cast it in the form q{x, t) = [J, Q{x, t)]. Since 
J is a regular element of f) then g(.x, t) E is determined 
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by |Z\| independent coefficient functions as follows: 



(9) 



where E±a are the root vectors of g and Z\+ is the set 
of positive roots of g, Z\ = Z\+ U (— Z\+). The root a is 
positive if a{J) > 0; by Ea, a E A and Hi, i = 1, . . . , r 
we denote the Cartan-Weyl basis of g with the standard 
commutation relations, see This choice of the gauge 
of I/(A) allows one to study the TV-wave type equations: 



i[J,Qt]-i[I,Q.] + [[I,Q],[J,Q]]^0- 



(10) 



Indeed eq. (|l^ ) allows Lax representation with the pair 
of Lax operators L{X) (H) and Af(A): 



M{X)i, = ^i- + [/, Q(x, t)] - XI j ^{x, t, X) = 0, (11) 

where / is also a constant regular element of / G i). 

There is a second canonical way to fix u p th e gauge of 
the Lax operator known as the pole gauge |13| : 



L^{x, t, X) = - XS{x, t)j xP{x, t, X) = 0, (12) 

where tpix, t, X) = g^^{x, t)ip{x, t, A), 



S{x,t) — KAg ■ J = g ^{x,t)Jg{x,t). 



(13) 



The gauge transformation which takes L{X) to L{X) 
g~^L{X)g{x, t) is performed with the Jost solution of ( 
taken at A = 0, i.e. g{x, t) E & and 

■ dg 



i— + q{x,t)g{x,t) = 0, 
dx 



lim g{x,t) = 1. (14) 
Applying the gauge transformation to M (A) we get: 

Mi>{x, t, A) = (^tj^ - Xf{S)^ Vi(x, t, A) = 0, (15) 

where f{S) is a function (in fact, a polynomial) to be de- 
termined below. Indeed only this choice of M (A) ensures 
the vanishing of the A^-term in the zero-curvature condi- 
tion [L, M] — 0; the terms proportional to A lead to the 
following form of the gauge equivalent NLEE: 



St 



dx 



f{S) - 0, 



(16) 



While the TV-wave type equations are well known their 
gauge equivalent ones to the best of our knowledge have 
not been derived yet. One of the difficulties in doing this 
is the necessity to express all factors in terms of S only. 

The gauge equivalent operators L{X) and L(A) have 
equivalent spectral properties and spectral data. This fact 
allows one to prove that the classes of NLEE related to 
L(A) and L{X) are also equivalent. 

In Section 2 we construct the NLEE gauge equivalent 
to the TV- wave equations ( |l0|) extending the results in Q . 



Namely we calculate the functions f{S) for the cases when 
g belongs to the classical series of simple Lie algebras. In 
Section 3 we briefly describe the interrelations between 
the scattering data of L(A) and L(A) and outline some 
of their reductions. We also reformulate the Riemann- 
Hilbert problem (RHP) for the gauge equivalent systems 
and describe the time evolution of the scattering data. In 
section 4 we extend the Zakharov-Shabat dressing method 
[0,||,|l^ for the gauge equivalent systems |^ and provide 
the general form of the 1-soliton solution of these system. 
These results are demonstrated on an example on the or- 
thogonal Lie algebra B2 ~ so(5) in section 5. This and 
other particular cases of Eqs. ( |l6[ ) describe isoparametric 
hypersurfaces |14|. 



2 General form of the gauge equivalent 
systems 

It is natural that f{S) = g^^(x, t)Ig(x, t), i.e., it is uniquely 
determined by /. Both J and / belong to the Cartan sub- 
algebra t) so they have common set of eigenspaces. 

The derivation of the corresponding functions f{S) is 
different for and for Br, Cr, series. Let first g ~ 
Ar = sl(n) with n = r + 1. Then 

J = diag(Ji, . . . , Jn), I = diag(/i, . . .,/„), 

and the only constraint on the eigenvalues Jk and Ik is 
trJ — trl — 0. The projectors on the common eigensub- 
spaces of J and / are given by: 

'rfc(J) = ITt—T 0,1,0,..., 0). (17) 

Next we note that 



/ = ^/fc7rfe(J). 



(18) 



fc=i 



In order to derive f{S) forg ~ sl{n) we need to apply the 
gauge transformation to (p^ with the result: 



(19) 



k=l 



i.e., f{S) is a polynomial of order n — 1. Obviously S is 
restricted by: 



Y[{S^Jk)=0, tr5^=trJ^ 



(20) 



fc=i 



for k — 2, . . . ,n. 

Let us now assume that g is an orthogonal or symplec- 
tic algebra. In the typical representation we can introduce 
the Cartan generators He^ which are dual to the orthogo- 
nal basic vectors in the root space. Each He^. has only 
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two non- vanishing eigenvalues equal to 1 and —1 respec- 
tively. Then we put: 

r r 



3 Fundamental analytic solutions (FAS) and 
scattering data for gauge equivalent systems 



k=l 



fe=l 



The direct scattering problem for the Lax operator 
based on the Jost solutions: 



Obviously the odd powers of H^^ also belong to g while 
the even powers do not. The projectors fk{J) onto H^^. 
then can be written down as: 



lini ipi^i 



iXJx 



1, 



lim 



(x, A)e 



iXJx 



Therefore 



fc=i 



(21) 



(22) 



and the scattering matrix: 

T{X) = {i,(x,\))-^cp{x,\). 



1,(28) 



(29) 



and applying the gauge transformation we get: 

r 

f{S) = g-\x,t)Ig{x,t)^Y.^kfk^^)- (23) 



fe=i 



Then the equation gauge e quiv alent to (^ is given by 
(|l|) with f{S) determined by (|2|). 

In addition 5(x, t) satisfies a set of nonlinear constraints; 
one of them is the characteristic equation: 



(24) 



fe=i 



where kq = if g ~ or Dr and kq = 1, if g ~ -Br- To 
construct the others we use the typical representation of 
Q. It this settings we see that all even powers of Hf.^ have 
trace equal to 2. Thus we have: 



The FAS ^^(x, A) of L{\) are analytic functions of A for 
A ^ and are related to the Jost solutions by ||l^ 

^^{x,\) = (t>{x,X)S^{X)^ip{x,X)T^{\)D^{X), (30) 

where T±(A), S^{\) and D^{\) are the factors of the 
Gauss decomposition of the scattering matrix: 

T(A) = T-{\)D+{\)S+{\) = T+iX)D-{X)S-{X). (31) 

Here S = S^^, the superscripts + (resp. — ) in T^{X) and 
S^{X) mean upper- (resp. lower-) triangularity. The diag- 
onal factors (A) are analytic functions of A for Im A > 
and Im A < respectively. 

On the real axis £,'^{x, A) and £,~{x, A) are related by 

e+(:E,A) =r(^,A)Go(A), Go(A) -^-(A)^+(A), (32) 

and the function Go (A) can be considered as a minimal 
set of scattering data in the case of absence of discrete 
eigenvalues of (||) jl5[0. 

If the potential q[x]t) of L{X) (H) satisfies equation 
( |l0|) then S^{X) and T^(A) satisfy the hnear equation: 



tr(j2'=) = 2^ =tr (5)2^ 



(25) 



''~dt 



X[I,S^ 







i — -- 

dt 



A[/,r=^]=0, (33) 



0=1 



for k — 1, . . . , r. The conditions ( p5| ) are precisely r inde- 
pendent algebraic constraints on S. Solving for them we 
conclude that the number of independent coefficients in S 
is equal to the number of roots \A\ of g. 

Both classes of NLEE possess hierarchies of Hamil- 
tonian structures. The phase space A4n-w of the N-wave 
equations is the linear space of off-diagonal matrices q{x,t); 
the hierarchy of symplectic structures is given by: 



while the functions D^{X) are time-independent. In other 
words D^{X) can be considered as the generating func- 
tions of the integrals of motion of (|l^) . 

In order to determine the scattering data for the gauge 
equivalent equations we need the FAS for these systems: 



^^ix,X)=g-'{x,m^{x,X)g- 



(34) 



(fc) 



i dxtr {SqAA''[J,Sq{xA)]) 



(26) 



The phase space A^s of their gauge equivalent equations 
( p^ ) is the nonlinear manifold of all S{x, t) satisfying equa- 
tions (p4), (pq). The family of compatible 2-forms is: 



where g- = linix ^ ^ oo g{x,t) = T{0). In order to ensure 
that the functions £,^{x, A) are analytic with respect to A 
the scattering matrix T(0) at A = must belong to the 
corresponding Cartan subgroup Sj. Then Equation ( |3^ ) 
provide the FAS of L. We calculate their asymptotics for 
X — > ±oo and establish the relations between the scatter- 
ing matrices of the two systems: 



= i I dxtr [SSAA" [S, SS{x, t)] 



lim ^+{x,X) =T{0)S+{X)T{0) 



(27) 



X — ' — oo 



— iXJxrp — 



Here A and A are the recursion operator of the TV-wave 
type equations (see [Q) and its gauge equivalent: A = 
g-^Ag{x,t). 



lim C{x,X) 

X — >oc 

with the result: 

f (A) = r(A)f (0) 



T-{X)D+{X)T{0) 



(35) 
(36) 

(37) 
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Obviously r(0) — 1. The factors in the corresponding If g ~ B,-.Dr the dressing factors take the form p6[ : 
Gauss decompositions are related by: 



u{x, A) = 1 + (ci(A) - l)Pi + (cr'(A) - l)P_i (45) 



S'±(A) = r(0)S'=^(A)f (0), f ±(A) = r±(A) 

7)±(A) = Z?±(A)f (0). (38) ^) = 1 + - 1 j ^1 + [^) - 1 j ^-^*6) 

On the real axis i+{x, A) and ^-{x, A) are related by: where P_i(x) = S'Pf (x)S'"\ Pi{x) is the rank 1 projector 

. . . . (H), P±i -g(-;P±ig(o)(x,t). Iffl-B, theniV-2r + l, 

e+(x,A)=r(a:,A)Go(A), Go{X) ^ S-{X)S+{Xp9) ^ W 

r 

with the normalization condition |(a;,0) = 1; again Go(A) 5" = ^(-l)''+^(Sfefe + E^^j.) + (-l)''£:,.+i^r+i; (47) 

can be considered as a minimal set of scattering data. k=i 

The numerous Z2-reductions that have been recently 

classified for the iV-wave equations ilfl , n7 | using the re- k — N ~ k + 1, {Ekm)ii = SikS„ii; if g ~ then N — 2r 
duction group introduced by Mikhailov |p_8||. They can eas- 
ily be reformulated for the gauge equivalent systems. Here 
we write down only two of them: 



S^J2(--^)'^'(^k-k+Ekk), (48) 



fc=i 



1) SHx,t) ^ KS{x,t)K-\ Ke^], K^ = l, 

J, If the dressing factors of the gauge equivalent equations 

2) S{x,t) — S{x,t) . (40) satisfy (44) then the projectors P±i satisfy the equations: 

Obviously each of the constraints 1) and 2) are compatible 

with equation ( [T6| ) and diminishes the number of indepen- i — - -|- Aj^_Pi5(o) — Aj^5(i)_Pi = 0, 

dent coefficients by a factor of 2. _ 

■ A+P_i5(o) - A+5(i)P_i = 0, (49) 

4 Dressing factors and 1-soliton solutions 

The main idea of the dressing method is starting from a A^ A^ d 

FAS |('^)(x, A) of L with potential 5(o) to construct a new S^i-) = 5(o) + i ^ ^ J --r{Pi{x) - P~i{x)). (50) 
singular solution ^^^^ (x, A) of the RHP (|39| ) with singu- 
larities located at prescribed positions Xf . Then the new The dressing factors can be written in the form: 
solutions ^^^^{x,X) will correspond to a potential of 

L with two discrete eigenvalues Xf. It is related to the u{x,X) — exp 
regular one by the dressing factors u{x, A): 



dx 

and the "dressed" potential can be obtained by: 



(51) 



if;Ax,X) = uix,X)it,.{x,X)uZHX), P(^) = ^1 - ^-1.^ ^^/i consequently S(x, A) be- 

'^'-i '^"1 longs to the correspondmg orthogonal group. 

u-{X) ^ lim u{x,X), (41) Making use of the explicit form of the projectors P±i(x) 

valid for the typical representation of B^. we have 

and the dressing factors for the gauge equivalent equations 
u{x, A) are related to u{x, A) by _ 2 



{x,X) ^ gj{x,t)u ^(a;, A = 0)u(a;, A)g(o). (42) 



{m\n) 



2 



^ (P„(x)i;„+P„„(x)£;_„), (52) 



If g ~ Ar then the gauge equivalent dressing factors are {m\n^ 

u(x, A) = 1 + f - 1^ Pi, ci(A) = ^Ll^ where we assumed 5(o) = J, g(o) = 1- Thus 
Vci(O) J A-Ai 

n{x)){m{x) ^^^^ hk{x,t) ^ no^kmo,ke^'''y'' -nQj,moj,e-^'''y'', 



^1^^) = / A ' (43) 
(TO(a;)|n(a;)) 



k=l 



where |no) and {mo\ are constant vectors and these dress- + ?^o,r+in^o,r+i, (53) 

ing factors satisfy the equation: / p f^,- ^ — o _ c 

I -^ks 7 ^ — 

du Pa = \ Pks, for a = e/c + Cs 

~ ^^{1)^ + ^^^{0) ^ 0- (44) iPfc,r+i, for a = efe 
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Here 1 < fc, s < r, /ii = Re Aj^, i^i = lmX]\ and 

yk = JkX + Ikt, Vk = -yk, yr+i = 0. (54) 

The corresponding result for the series is obtained for- 
mally if in the above expressions (|5^) , (|5^) and we put 

no.r+i = TOo.r+i = 0. Thus Pk.r+i = Pr+1 fc = and the 
last term in the right hand side of {mln) (|5^) is missing. 

The reductions ( |4C| ) apphed to the 1-soliton solution 
constraint the vectors~lno), (too| and the eigenvalues : 

1) \no)=K\m*), Ar = (A+r, 

2) \no} = |mo), A^ = -A+. (55) 

The iV-soliton solutions can be obtained by applying suc- 
cessively N times the dressing procedure. 

5 Example: g ~ B2 algebra 

This algebra has four positive roots: ei ± 62, ei and 62. 
The corresponding 4-wave system has the form: 

i(Ji - J2)qio,t - - h)qw,x + 2Kq-i_iqQ^ = 0, 

iJiqoi.t ~ ihqoi,x + K{qnqi2 + qnqlo) = 0, 
iJiqii.t - ihqii.x + K{qi2qoi - qwqoi) = 0, (56) 

i( Ji -I- J2)qi2.t - + h)qi2,x - 2Kqii(7oi = 0. 

where k = J1/2 — J2-fi and the subscripts 10, 01, 11 and 
12 refer to the roots 61 — 62, ei, 62 and 61 -I- 62 respectively. 
This system has applications in nonlinear optics [^0] and 
in differential geometry [|l^. Its gauge equivalent is: 

St - hSx - h{s% = 0, 



/i 



I2J1 — I1J2 



/3 = 



I1J2 



JlJ2(Ji — J2) JlJ2{Ji — J2) 

where the 5x5 matrix S is constrained by: 

trS^ = 2(j2 + J2), tr5^ = 2(jf -I- J^), 

5(52- J2)(52- J|)=G. (58) 

We write down the 1-soliton solution for a special choice 

"0,1 = 1, "-0,2 = P, no^3 = VAp'^ - 1) = 
no,k = "o,fc' "^o,fc = no^k- (59) 

of the soliton parameters with p > 1 and real. The choice 
( Isoj ) satisfies dssl) w ith K — 1. Inserting it into the general 
formulae (^-(^4|) we get P_q = P* with: 

(m|n) = 2(sinh^ i/iyi + cosh^ '^iy2), 

hi ~ sinh2i^i?;i, ft,2 ~ sinh 2i/i?;2, 
P,^±,^ = pe-*^i(^i±^^)coshi.i(yi Ty2), (60) 

k 



P. 



v/2(p2 _ i)e-wy. sinhi.i2/fc. 



1,2. 



If p = 1 we get a 1-soliton solution associated with the 
D2 ~ © subalgebra; p = gives a 1-soliton solution 
for the so(3) subalgebra of ^2. In both subcases the sub- 
sets of roots ({±ei ± 62} and {±ei, ±62} resp.) for which 
Pa ^ contain only roots with the same length. 



6 Discussion 

We derived the explicit form of the NLEE gauge equivalent 
to the iV-wave equations related to the classical simple Lie 
algebras g. These equations are Hamiltonian ones. Their 
phase space is a nonlinear one, since S is restricted by 
the nonlinear constraints (|2j), (pq). 

The gauge covariant formulation of the spectral de- 
compositions of the recursion operators A and its gauge 
equivalent A allows one to expect that the hierarchies 
of symplectic structures will satisfy in analogy with the 



NLSE-HFE case the relation [2 



(fe) 



[2 



(fc+2) 



'N-w — 

Another open problem is the study of the Z2-reductions 



of (|16|) along the ideas outlined in 
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